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Abstract 

This paper explores the consequences of non-Gaussian cosmological perturbations for the forma- 
tion of primordial black holes (PBHs). A non-Gaussian probability distribution function (PDF) of 
curvature perturbations is presented with an explicit contribution from the three-point correlation 
function to linear order. The consequences of this non-Gaussian PDF for the large perturbations 
that form PBHs are then studied. Using the observational limits for the non-Gaussian parameter 
/jvx, riew bounds to the mean amplitude of curvature perturbations are derived in the range of 
scales relevant for PBH formation. 
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I. INTRODUCTION 



The non-linear theory of perturbations is an important part of the study of the early uni- 
verse because only at this level of precision one can distinguish specific models of inflation 
and even rule out some of them through observations of the correspondent non-Gaussian 
statistics. The possibility of observing non-Gaussian correlations in the temperature fluctu- 
ations of the cosmic microwave background radiation (CMB) [l| and in the count of large 



scale structure (LSS) has fostered the study of 

models of the origin of structure (see e.g. Refs. 
review of both observational and theoretical aspects of non-Gaussianity,). 

The properties of the early universe imprinted in the observed perturbations are better 
studied when we make use of the curvature perturbation TZ(t, x) on comoving hypersurfaces 
[55] . This quantity is defined in a gauge- invariant way and, when sourced by adiabatic matter 
perturbations, the growing perturbation mode remains constant on scales larger than the 
cosmological horizon H~ l , where H = dln(a)/dt is the Hubble parameter P. fiol] . 

A fundamental observable of cosmological perturbations is the mean amplitude. The 
classical amplitude of metric perturbations 1Z C \ is derived by solving the perturbed Einstein 



ay oi nc 

0,y,3 



non-linear perturbations for distinct 



and see Ref. 



8] for an extended 



equations to linear order. Statistically this classical amplitude is written in terms of the two 
point correlation function as 

(^ G (k 1 )^ G (k 2 )) = (2 7 r) 3 5(k 1 + k 2 )\1l cl (k)\ 2 , (1) 

where Tlcfe) are Gaussian perturbations in Fourier space, and ki and k 2 are vectorial 
Fourier momenta. The two-point correlator defines also the dimensionless power-spectrum 



[55] In the longitudinal gauge, the scalar isotropic curvature perturbation or Bardeen potential $ modifies the 
Friedmann background metric to the form: 

ds 2 = - [1 + 2$(t, x)] dt 2 + a(t) 2 [1 - 2$(i, x)] dx 2 . 

where t is the cosmic time, a(t) is the scale factor and <ix is the comoving spatial part. The comoving 
curvature perturbation TZ{t, x) is defined in the comoving gauge from the metric 

ds 2 = -N 2 (t,x)dt 2 + a(£) 2 e 2K(t ' x) 5y (dx* + N i (t,x)dt) (dx j + N j (t,x)dt) , 

where dij is the Kronecker delta, and N and iV l are functions algebraically determined from the symmetries 
of the perturbed equations of motion. 
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V(k) as 



2tt 2 

(^ G (k 1 )^ G (k 2 )) = (2vr) 3 5(k 1 + \i 2 )—V{k x ). (2) 



The power-spectrum encodes important information of the underlying cosmological 
model. For example, for perturbations generated from a single inflationary field <fi with 
a potential V dominating the cosmological dynamics, TZc{k) is a random field of perturba- 



tions about the quasi-de Sitter background, and the power-spectrum is given by |12|] 



H 4 V 3 

where mpi = 8nG is the Planck mass and where an asterisk denotes evaluation at t*, the time 
when the relevant perturbation mode exits the cosmological horizon, k = a^H* = a(t*)if(t*). 

The tilt of the power-spectrum is parametrised with a second observable, the spectral 
index which is defined as 

n, = ^V(k). (4, 

With this definition, n s < means that the power-spectrum is larger for large scales, in such 
case we have a red spectrum. Equivalently n s > corresponds to a larger power for small 
scales and this is called a blue spectrum. 

The power-spectrum and the tilt are derived directly from linear perturbations. Assuming 
linearity for the perturbations of the CMB, it is possible to determine with great accuracy 
the numerical values of the power-spectrum and its tilt on scales larger than the horizon at 
the time of last scattering. At such scales, the most recent observations of the CMB and the 
galaxy count give V = 2.4 x 10~ 9 , n s = —0.05 ± 0.01 [1]. Any successful theory of structure 
formation must meet these values at the relevant scales. Higher order correlations of the 
perturbation field 1Z offer an exciting way to distinguish between these cosmological models. 

The deviations from Gaussianity are described to lowest order by the non-linear parameter 
/nl- Mathematically, this parameter appears in the expansion 13] 

n{k) = n G (k) - ^hL(K G *n G (k) - (n 2 G )), (5) 

where a star denotes a convolution of two copies of the field. We use this definition through- 
out the paper. At the present time, the precision of the observations can only provide limits 
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for /ml. The WMAP satellite gives the constraints —54 < /nl < 114 at the 95% confidence 



level 



!]• 



When one computes second-order perturbations to t 



re Einstein equations, higher order 



correlations become non-vanishing (see e.g. Refs. QQ). Alternatively such correlator 
can be obtained from the third-order quantum perturbations to the Einstein-Hilbert action 
(Pioneering works using this method are Refs. H Q,Q). The three-point correlation 
function in Fourier space at the tree level can be directly computed from the definitions in 
Eqs. CO)© and © to find Q 
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(n{k y )n{k 2 )n{h)) = - (2tt) 3 <5 4vr 4 -/ NL 



V{ki)V{k 2 ) 



+ 2 perm . 



Non-linear correlations have been explored not 



also in counts of galaxies and clusters 



20 



only in the CMB anisotropies 



HQ, 



(6) 



but 



211 ] . The latter deserve special mention. 



The distribution of galaxies and clusters account for density perturbations that formed 
such structures in the matter-dominated era. The statistics of these objects gives us extra 
information about the primordial non-Gaussianities in the universe. However, in studying 
the formation of galaxies and clusters from the seed perturbations, one must be careful in 
the treatment of intermediate stages because the growth of perturbations after horizon entry 
is intrinsically non-linear. This may blur the primordial non-Gaussian statistical effects (For 
a review, see Ref. jsj). 

In this paper we study non-Gaussianity for a range of wavelengths not covered by the ob- 
servations mentioned above by looking at the formation of Primordial Black Holes (PBHs). 
PBHs were formed during the radiation era of our universe from large amplitude pertur- 
bations. The importance of non-Gaussianity in this context is that a small deviation from 



22 



he effects of non-Gaussian 



23 



24J but a precise quan- 



linearity would change the probability of formation significant 
perturbations on PBHs have been studied for specific models 
tification of the non-Gaussian effects is still required and it is only now, with a much better 
understanding of the effects of higher order perturbations, that we are able to describe the 
general effects on PBHs. The discussion of non-Gaussian effects on PBH formation turns 
crucial in the light of recent works (e.g. 25J) which claim that only through exotic exten- 
sions of the canonical slow-roll inflationary potentials can produce considerable amounts of 
PBHs. Here we explore whether non-Gaussian perturbations could ease this dificulty. 

To look at the effects of non-Gaussianity we make use of a non-Gaussian PDF for cur- 



vature perturbations with a direct dependence on the non-linear parameter /nl as derived 



in 



26[. We adapt this PDF to calculate probabilities for pertrbations relevant to PBH for- 



mation. We also explore the amount of non-Gaussianity that is allowed from the bounds to 
the total mass fraction of PBHs. In principle, the PBH abundance could be used to probe 
non-Gaussianity on scales not reached through CMB or LSS observations. We find that even 
though the possible limits on /nl in the relevant mass regime are weak, its current allowed 
values in the CMB scale can be used to generate new bounds on the mean amplitude of 
curvature perturbations 1Z for the wavelength range < X/Mpc < 1. 

This paper is organised as follows. In Section [III we construct the appropriate non- 
Gaussian probability distribution function (PDF) to calculate the probability of non-linear 
perturbations forming PBHs. In Section II III we characterise the effects of the non-Gaussian 
PDF on the probability of PBH formation. We recover the results presented in previous 



studies 



22 



23| and reconcile the discrepancy in their conclusions. In Section HVl we make use 



of the Press-Schechter formalism to derive the abundance of PBHs including non-Gaussian 
effects. We show how the constraint on the amplitude of perturbations is modified by the 
introduction of a non-Gaussian mass function and report on new bounds on the variance of 
curvature perturbations. We conclude in Section [V] with a summary of results. 



II. THE NON-GAUSSIAN PDF 



In linear perturbation theory one makes use of the central limit theorem to construct 
the probability distribution function (PDF). To first order, the perturbation modes are 
independent of each other. If we define the field of perturbations TZ with zero spatial 
average, then the central limit theorem indicates that the PDF of TZ is a normal distribution 
dependent on a single parameter, the variance (TZ 2 ), 

An immediate effect of including the non-linearity of the perturbations is the interaction 
of distinct perturbation modes. Indeed, as shown in Eq. (jSJ), the parametrisation of higher 
order contributions to TZ is described in terms of convolutions in Fourier space. The require- 
ments of the central limit theorem are not met, so it can no longer be used to construct the 
PDF. 
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In a recent paper a first-order correction to the Gaussian PDF was explored, and a new 
PDF was derived, which includes the linear order contribution from the 3-point function. 
In the following we describe the elements of such derivation. 

As is customary in the treatment of perturbations, we first smooth the field over a given 
mass scale ku with a window function Wu{k), 

dk 



(2tt) 



(8) 



Here we use a truncated Gaussian window function: 



W M (fc) = ©(Aw - fc) exp ( - u , 



k 2 
k 



M 



(9) 



where is the Heaviside function and the fiducial scale /c max is introduced to avoid ultraviolet 
divergences. The amplitude of the perturbation is parametrised by the central value of the 
configuration 7£(x = 0) = C,n. This parametrisation is particularly useful when we pick 
the relevant perturbations for the formation of PBHs 27J. The non-Gaussian probability 
distribution function for a perturbation with central amplitude is 26] 



Y 3 



Y 3 



(10) 



where Y\ is the variance of the smoothed field, related to the power-spectrum by 



Y\{M) 



fw^k)V(k), 



(11) 



and the factor J encodes the non-Gaussian contribution to the PDF: 



dki dk 2 dk 2 



(2tt) 9 
dki dk 2 dk 3 



W M {k 1 )WM{k2)W M {h){n{k 1 )'R{k 2 )'R{k z )) 

V{k x )T{k 2 \ 



-AY,*) ^ 



hk 2 



2 perm. 



(12) 
(13) 



where the last equation is valid at tree level in the expansion of (1Z1Z1Z). This limitation is 
justified as long as the loop contributions to the three-point function, generated from the 
convolution of 7?.-modes, are sub-dominant. Such is the case for /nl ^ f° r ^ values 

of k. As will be seen in section HV] the tested /nl does not exceed such values. 

Further details of the derivation of the PDF in Eq. ( JTOl) can be found in Ref. 26]. 
Here is sufficient to say that the time-dependence of this probability is eliminated when 
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the averaging scale is kyt < a(t)H(t) because the growing mode of the perturbation 1Z is 
constant on superhorizon scales jijl. [lo|. 

We now proceed to integrate the non-Gaussian factor in Eq. (j!3p over the amplitudes 
and scales relevant to PBHs. The expression in Eq. ([TBI) is integrated within the limits /c m ; n 
and /c max defined conveniently to cover the relevant perturbation modes for PBH formation. 
These objects are formed long before the matter-radiation equality, so in the large-box (small 
wavenumber) limit of the integral ffl3|) we can choose k m \ n = Hq, the Hubble horizon today. 
This is a reasonable lower limit for integrating perturbations relevant for PBH formation 
. At the other end of the spectrum, the smallest PBHs have the size of the Hubble horizon 



at the end of inflation. We therefore use k max = a, (£end ) -^end ? the comoving horizon at the 
end of inflation, as a suitable upper limit. It is important to mention that, even though the 
integral in Eq. ffl3|) should add all /c-modes, finite limits are imposed to avoid logarithmic 
divergences. Moreover, due to the window function Wu(k), the dominant part of the integral 
is independent of the choice of integration limits as long as they remain finite. 

We solve integral (fTBl for the limit of equilateral configurations of the three-point cor- 
relator, that is, correlations for which k± = k% = k%. This is not merely a computational 
simplification. In the integral, each perturbation mode has a filter factor Wyi{k) which, 
upon integration, picks dominant contributions from the smoothing scale kyi common to 
all perturbation modes [56]. We take this argument as an ansatz, in which case J can be 
written in the suggestive way: 

3 ' = -\ C 't^»« p <*>i'(M ' (14) 

In the following sections we compute J numerically for inflationary perturbations gener- 
ated in a single field slow- roll inflationary epoch, and for the case of constant /nl- This will 
be used to test the effects of non-Gaussianity on PBHs. 



[56] This is an important restriction to the kind of non-Gaussianities that can be tested by PBHs [29(. If 
a strong non-Gaussian signature is encountered exclusively for isosceles triangulations, this cannot affect 
the formation of PBHs. 
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III. NON-GAUSSIAN MODIFICATIONS TO THE PROBABILITY OF PBH FOR- 
MATION 



The simplest models of structure formation within the inflationary paradigm are those 
where a single scalar field drives an accelerated expansion of the spacetime and its quantum 
fluctuations give birth to the observed structure in subsequent stages of the universe (For an 
extended review of the inflationary paradigm, see Ref. 30[). Here we study the curvature 
perturbations generated at the time of single scalar-field inflation including contributions 
from non-linear perturbations. 

The effects of non-Gaussianity on PBHs have been explored in the past but with incon- 
clusive results. Bullock and Primack 2^] studied the formation of PBHs numerically for 
perturbations with blue spectra (n s > 0) and non-Gaussian contributions. The motivation 
for this was that any inflationary model with a constant tilt and consistent with the nor- 
malisation of perturbations at the CMB scale, must have a blue spectrum to produce a 



significant number of PBHs 31 



32l |. Their analysis is based on the stochastic generation of 



perturbations on superhorizon scales, together with a Langevin equation for computing the 
PDF. For all the cases tested, the non-Gaussian PDF is skewed towards small fluctuations, 
so PBH production, which integrates the high amplitude tail, is suppressed with respect to 
the Gaussian case. An example of the kind of potentials studied in Ref. 22] is 



Vi(<f>) 



1 + arctan 



( -*-) , for 



1 + ( 4X1 ° 33 ) (4) ' fOT 



>0, 
< 0. 



(15) 



where Vo is the amplitude of the potential at = 0. 

Another way of generating large perturbations in the inflationary scenario is to consider 
localised features in the potential dominating the dynamics. As one can see from Eq. (j3J), 
an abrupt change in the potential would generate a spike in the spectrum of perturbations. 
The effects of non-Gaussianity for an inflationar y m odel producing features in an otherwise 
red spectrum (n s < 0) were explored by Ivanov 23| using the toy model 



W) 



A- 



for < 0i, 



< 



A(0 2 - 0) + A% for 2 > > l5 



(16) 



for > 2 . 

where A and A are coupling constants. Through a stochastic computation of the PDF, it 



S 



was found that large amplitude perturbations were more abundant for a non-Gaussian PDF 
than for a Gaussian one. 

To understand this difference and generalise the effects of non-Gaussianity, we look at 
the fractional difference of the Gaussian and non-Gaussian PDFs: 



P 



NG 



P, 



G 



Pf 



r 3 



n 



J_ 

Y 3 



(17) 



Note that both Refs. 



22] and |23j use perturbations generated in a piecewise slow-roll 



inflationary potential for which inflation is controlled by keeping the slow-roll parameters 
e = l/2(mpiV /V) 2 and rj = mp X {V" /V) smaller than one. Here we make use of the slow- 
roll approximation to explore the qualitative effects of Eq. (fTTj) . To linear order, there is a 
straightforward expression for the spectral index in terms of these parameters 121 ]. 



n s = 2(77 - 3e). 



(18) 



On the other hand, carrying a first order expansion on slow-roll parameters, Maldacena, 
jrovides an expression for the non-linear factor /nl in terms of the slow-roll parameters 

m 



/nl =4 K + Hk)n t ) = I (77 - 3e + 2F{k)e) 
12 



(19) 



where n t = 2e is the scalar-tensor perturbation tilt and T{k) is a number depending on the 
triangulation used. For the case of equilateral configurations, when T = 5/6, 



f 5/4 
/nl =-Jv- r 



(20) 



eq 



We use this last result to evaluate the integral (fT4l) . The non-Gaussian effect on the PDF 
is illustrated in FigJT] for the potentials given by Eqs. ffT5l) and ffTBT) in terms of the frac- 
tional difference (TTTj) . This last factor represents the skewness of the non-Gaussian PDF. 
Consequently the sign of /nl is what determines the enhancement or suppression of the 
probability for large amplitudes in our non-Gaussian PDF. Indeed, the non-Gaussian 
contribution encoded in the factor J is the sum of the /nl value over all scales relevant for 
PBH formation. For the two cases illustrated, the scalar tilt n s dominates over the tensor 
tilt n t in a way that the sign of /nl incidentally coincides with that of n s . 



9 



IV. CONSTRAINTS ON NON-GAUSSIAN PERTURBATIONS IN THE PBH 
RANGE 



Primordial Black Holes are objects that collapsed from large-amplitude perturbations 



at times previous to photon decoupling 



33j. The energy density of PBHs formed during 



inflation is diluted by the superluminous expansion. In consequence, a significant production 
of PBHs can only take place after inflation. Depending on the model of inflation, PBHs can 
cover a wide range of masses Mend ~ 10~ 48 M Q < Mpbh < 1OM at scales much smaller 
than galaxies M ~ 10 9 M Q or clusters M ~ 10 13 M Q . We will call this the PBH mass range. 
The lack of direct observations of PBHs limits their cosmological abundance. Assuming 
they are nonetheless present, there are three ways in which PBHs can affect the evolution of 
our universe and converselly, we can impose constraints to the abundance of such objects. 
First, the current density of PBHs cannot exceed the amount of dark matter density, i.e., 
Qpbh(M > 10 15 g) < tt DM = 0.47 (with the WMAP-III value taken at 3a confidence level 
ll). Second, the Hawking radiation 34J from PBHs, can also generate the radiation observed 




0.2 



0.6 



c 



FIG. 1: The fractional departure from Gaussianity is plotted for two types of non-Gaussian distri- 
butions Png (fTO]) . For the potential in Eq. ([15]) , / NL > 0. The potential of Eq. (Tl6|) gives /nl < 
and its correspondent PDF is shown by a dashed line. 
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at various wavelengths in our universe 35[. The lifetime t evSiP for PBHs is a function of its 

mass 



t 



cva] ) 



1.2 x 10 



M 
mpi 



sec. 



(21) 



From this relation we immediately infer that PBHs of mass M cvap = 5 x 10 g are evaporating 
today and the observation of the gamma-ray background constrains their present density 



parameter to f2pBH(^eva P ) ^ 5 x 10' 



35 



36 



37 



This is the tightest constraint on the 



abundance of PBHs. A third cathegory of constraints is relevant for lighter PBHs. Black 
holes with mass M < M evap have already evaporated and the decay products should not 



spoil the well understood chemical history of our universe (see e.g. Refs. |39|, |40j). 

To calculate the PBH mass fraction we make use of a standard Press-Schechter formalism 



4=11 ] . This formula integrates the probability of PBH formation over the relevant matter 
perturbation amplitudes 33]]. This is interpreted as the mass fraction of such PBHs at the 
time of formation, 



/9pbh(> M) 



= 2 / ¥(5(M))d5(M) 



th 



-cxp 



si 



2af(M) 



(22) 
(23) 



Here 5 = 5p/p is the matter density perturbation erf is the corresponding variance and 5th is 
the threshold amplitude of the perturbation necessary to form a PBH. By integrating over 
a smoothed perturbation, this integral is equivalent to the mass fraction of PBHs of mass 
M > ^' 2 M H « 7 3 / 2 fc M /(27r) [33j, where 7 is the sound-speed squared at the time formation. 

is valid only for a Gaussian PDF. 



Note that the approximation 

The integral (I22p establishes a direct relation between the mass fraction of PBHs and the 
variance of perturbations. The set of observational constraints on the abundance of PBHs is 



listed in Tab 
cosmologies 



31 



and has been used to bound the mean amplitude of 5 defined by for distinct 



32 



42 



431 ] . The Press-Schechter formula has also been tested against other 



methods such as peaks theory 

The threshold value 5$, used in Eg . (|22l) has changed with the improvement of gravita- 



tional collapse studies 



27 



33 



45 



46 



471 ] . Here we use the value 5 t h = 0.3 for convenience 



[57]. The corresponding threshold value of the curvature perturbation can be deduced from 



[57] Recent studies suggest that this value is dependent on the profile of the curvature perturbation 
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the relation [3 



S k (t) 



2(1 + 7) 



k 



(24) 



5 + 3 7 

which at horizon crossing during the radiation dominated era gives IZth = 0.7. Note that 
the high amplitude of the perturbations relevant to PBH formation necessarily require a 
non-linear treatment of their statistics. This is the major motivation for our analysis. 

We adopt the Press-Schechter formula derive the non-Gaussian abundance of PBHs. The 
use of the Press-Schechter integral on distributions of curvature perturbations is not new. 
In Ref. j|3(| the Press-Schechter formula is used to integrate curvature perturbations which 
never exit the cosmological horizon. We apply the integral formula in Eq. (1221) to the non- 
Gaussian probability distribution (fT0|) . The result of the integral is the sum of incomplete 
Gamma functions r; nc and an exponential: 



1 



r in c f 1/2, 
J J_ 

v^^(M) 



2 T- 



c 



2Y? n {M) 



3 exp 



(25) 



The Taylor series expansion of these functions around the limit E^/Cft = gives 



P(M) 



CthV^ 



cxp 



2S^(M) 



1 



Cth 



J_ 

Y 3 



Cth 



1 



(26) 



For the mass fraction shown in Eq. 02 6p . the observational limits of Table I could in 
principle constrain the values of the variance and of /nl- However, when we normalise 
the amplitude of perturbations to the observations at CMB scales, the limiting values for 
/nl are too large, of order /nl ~ 10 4 , and thus inconsistent with perturbation theory at the 
level of the expansions performed in this paper. In fact, the expansion in Eq. (jSJ) shows that 
when 



|/nl| > ^ 



1 



V{k) 



203.2, 



(27) 



the cuadratic term of Eq. (jHJ) dominates over the linear term, and in the computation of the 
three-point function Eq. ([6]), the loop contributions become dominant. The computation 



Furthermore some early universe models may not generate the profiles required for PBH formation. This 
is a crucial aspect of PBH formation currently under investigation [49 1 
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of non-Gaussianities in this case goes beyond the scope of this paper and should be treated 



carefu 
Refs. 



ly elsewhere, (for discussions on the origin and magnitude of the loop corrections see 



51 



52 



sal) 



Here we restrict ourselves to use the marginal values allowed for /nl from WMAP-III ob- 
servations and look at the modifications that large non-Gaussianities bring to the amplitude 
of perturbations at the PBH scale. 



Or 




FIG. 2: The constrains in Table I are plotted together with only the smallest value considered for 
every mass. 



In Fig. [2] we plot the set of bounds to the initial mass fraction of PBHs listed in Table I. 
The corresponding bounds on are shown in Fig. [3] for the Gaussian and non-Gaussian 
cases. Independently of the model of cosmological perturbations adopted, one can use 
the observationa limits of /nl to modify the bounds for on small wavelengths. The 
tightest constraints on comes from perturbations of initial mass M ~ 10 15 g. With the 
non-Gaussian modification the limit is log < —1.2, compared to the Gaussian case 
log (S^) < —1.15. As shown in Fig. [3l the modification to cannot be much larger if we 
use the limit value of Eq. ( 1271) . 
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V. FINAL REMARKS 



In this paper we have computed, to lowest order of non-linearity, the effects of non- 
Gaussian perturbations on PBHs formation. 

We use curvature perturbations with a non-vanishing three-point correlation to find an 
explicit form of the non-Gaussian PDF with a direct contribution from the non-linear pa- 
rameter /nl- We have shown how the sign of this parameter determines the enhancement 
or suppression of probability for large-amplitude perturbations. Using the simple slow-roll 
expression for /nl in the context of single field inflation, We have resolved previous discrep- 
ancies in the literature regarding effects of non-Gaussianity on the abundance of PBHs. 

As a second application of the non-Gaussian PDF we have employed the Press-Schechter 
formalism of structure formation to determine the non-Gaussian effects on PBH abundance. 
We have shown how the PBH constraints on the amplitude of perturbations can be modified 
when a non-Gaussian distribution is considered. The maximum variance allowed by 



-1 




FIG. 3: A subset of the constraints on T-ji from overproduction of PBHs is plotted for a Gaussian 
(black line) and non-Gaussian correspondence between (3 and T-ji, equations (|23p and (|26p respec- 
tively. The green dashed line assumes a constant /nl = —54 and the blue dotted line a value 
/nl = 
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PBH constraints is sensitive to /nl, producing the limit E^(M = 10 15 g) < 6.3 x 10~ 2 for 
/nl = —54. This limit is, however, much larger than the observed amplitude at CMB scales, 
where « 4.8 x 10~ 5 . The order of magnitude gap between the mean amplitude observed 
in cosmological scales and that required for significant PBH formation remains almost intact 
and, as a consequence, non-Gaussian perturbations do not modify significantly the standard 
picture of formation of PBHs. 
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Table I Constraints on the fraction of PBH density to the density of the universe as described in 



Ref. 
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Qpbh (today) < 1 
X-rays from 

evaporating PBHs 

evaporated PBHs X-rays 

pair production at nucleosynthesis 
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